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GLOBAL  PNS  SOLUTIONS  KOR  LAMINAR  AND  TliRBUl.KNT  BLOW 


jV  ifk 

S.C.  Rubin  and  D.R.  Reddy 

University  of  Cincinnati 
Cincinnati,  Ohio  A 5221 


Abstract 

A  multi-sweep  relaxation  procedure  is  applied 
for  inviscid  and  parabolized  (pressure-elliptic) 

Nav ier-Stokes  (PNS)  equations.  Boattail,  finite 
flat  plate  and  NACA  0012  airfoil  geometries  are 
considered  for  incompressible  and  subsonic  inviscid, 
laminar  and  turbulent  flow.  The  equations  are 
written  in  a  conformal  body  fitted  coordinate  frame 
and  differenced  on  a  staggered  grid  in  order  to 
give  second-order  accuracy  for  the  inviscid  flow 
and  somewhere  between  first  and  second-order 
accuracy  for  the  PNS  solutions.  A  full  second- 
order  scheme  is  also  discussed.  Separation, 
trailing  edge  and  stagnation  point  flow  are  evalu¬ 
ated.  The  effects  of  normal  pressure  gradients 
for  laminar  and  turbulent  flows  are  compared.  A 
multi-grid  procedure  is  applied  in  order  to  speed 
convergence  rates  for  fine  meshes  and/or  large 
computnt ional  domains. 

1 .  Int  roduc  t  ion 

The  reduced  form  of  the  Nav ier-Stokes  equations 
originally  termed  parabolized  Navier-Stokes  (PNS) 
and  more  recently  semi-elliptic,  partially  parabolic 
or  thin  layer  are  considered  herein  bv  relaxation 
methods  for  the  evaluation  of  incompressible  and 
subsonic  flows  with  strong  pressure  interaction 
and/or  separation.  Previous  analysis  bv  the 
authors  and  solutions  for  laminar  incompress ible 
flow  for  trough  and  boattail  configurations  are 
given  in  references  1-4,  for  subsonic  and  transonic 
flow  over  a  cone-cyl inder-boat ta il  geometry  in 
references  6,  7  and  for  supersonic  flow  over  a 
cone  in  reference  8.  The  present  procedure  is 
applicable  to  inviscid  as  well  as  viscous  inter¬ 
acting  flow  and  can  be  classified  as  intermediary 
between  interacting  boundary  layer  theory  and  con¬ 
ventional  Navier-Stokes  or  PNS  technique. 

The  PNS  equations  were  first  applied  for  hyper¬ 
sonic  problems'  where  the  contribution  of  the 
pressure  gradient  p  in  the  longitudinal  (x)  momen¬ 
tum  equation  is  negligible  and  can  be  neglected. 

The  svstem  is  then  mathematically  parabolic  and 
can  be  solved  as  an  initial  value  problem  by  march¬ 
ing  techniques.  For  lower  Mach  numbers,  where  px 
must  be  retained,  an  elliptic  influence  associated 
with  the  pressure  interaction  through  the  subsonic 
portion  of  boundary  layer  appears  in  the  PNS 
system. 10,11  Single  sweep  marching  then  leads  to 
an  ill-posed  initial  value  prob  lem* **  and  expo¬ 

nentially  growing  departure  solutions  appear  for 

step  sizes  Ax  (Ax)  .  ,  where  (Ax)  .  is  propor- 
r  min’  min 

t ional  to  the  extent  of  the  subsonic  portion  of 

the  flow.®* 


*  Head  and  Professor,  Dept,  of  Aerospace  Engineer¬ 
ing  and  Appl  it'd  Mechanics,  Member  AIAA. 

**Kese«irch  Associate,  Dept,  of  Aerospace  Engineer¬ 
ing  and  Applied  Mechanics. 


Alternatively,  an  eigenvalue  analysis  has 
shown  that  the  px  term  can  be  structured  so  that 

P  ~  .i(p  )  +  (l-.)(p  )  ,  where  .  (M)  is  a  function 

x  x  p  x  e 

of  a  local  Mach  number  M  and  0  .  .  1 ,  -(0)  =  0, 

.(1)  -  1.  The  term  . (p^)  represents  the  parabolic 

(marching)  contribution  to  p^_  and  (!-«.)  (p^)  repre¬ 
sents  the  elliptic  contribution  to  p  .  This  latter 
term  contains  the  upstream  influence.  If  included 
in  a  marching  calculation  where  the  M  1,  the  ill- 
posedness  of  the  initial  value  problem  is  once 
age  in  apparent  with  the  appearance  of  the  departure 
solutions  that  reflect  different  downstream 

boundary  conditions  for  the  pressure.  The  (Ax)  . 

m  in 

and  ..(p  )  descriptions  can  he  reconciled  by  defin- 
x  P  _ 

ing  Ax  as  Ax/,  in  the  general  marching  analysis. 

For  given  Ax,  o:(M)  and  therefore  Ax  vary  across 
the  subsonic  layer  in  order  to  obtain  departure 
free  marching  solutions.  This  procedure  is  incon¬ 
sistent  for  flows  where  is  important. 

In  recent  papers  by  the  authors,^  ^  it  has 

been  shown  that  for  incompressible  flow,  (Ax) 

.  ,  ^  min 

is  proportional  to  the  total  extent  of  t  hi* 

computational  boundary  y«  in  the  surface  normal 

direction.  This  implies  that  for  vkl  »  (Ax)  . 

M  min 

or  with  .(0)  =  0  the  entire  px  contribution  is 
elliptic  and  therefore  must  be  neglected  for  depar¬ 
ture  free  single  sweep  marching  procedures.  For 
subsonic  flows,  o  •  1  over  a  range  of  y  values  so 

that  (Ax)  >>  1.  However,  from  the  p  split  of 
min  *  *  x  ‘ 

reference  13,  shown  previously,  it  is  seen  that 

(Ax)  .  for  subsonic  flow  is  reduced  bv  the  factor 
min 

(1-y)  over  the  value  obtained  in  references  2-5  for 

incompress ib] e  flow,  i.e.,  (.’x)  —  (l-~)(..x)  • 

’  *  min  min/M=0 

In  order  to  circumvent  the  i 1 1 -posedness  of 
single  sweep  PNS  methods,  a  gl oba 1  pressure  relaxa- 
Jt^ion  or  multiple  marching  procedure  has  been  pro¬ 
posed  bv  the  authors^*"®  for  the  entire  Mach  number 
r an?at> »  i»e.,  incompressible  t o  sup e rso n  i c  f l_ow . 

This  requires  an  appropriate  and  exact  "forward" 
or  mid-point  difference  treatment  of  the  p^  con¬ 
tribution  since  this  term  alone  (for  attached  flow) 
contains  the  upstream  or  elliptic  influence. 

S  ign i  f  ic an 1 1  y ,  onl  y  the  _p re s k  u re  (andj/e  1< >r  i  t  j_es  in 
in  regions  of  reversed  flow)  must  be  stored  during 
the  relaxation  process.  This  leads  to  a  signifi¬ 
cant  reduction  over  conventional  Navier-Stokes  and 
PNS  methods. 


Consistent  (Ax  arbitrary),  departure  free 
(Ax  -►  0)  and  rapidly  convergent  solutions  have 
previously  been  obtained  for  laminar  incompressi¬ 
ble, subsonic,^  transonic'7  and  weaklv  inter- 
8 

acting  supersonic  flows.  Strong  pressure  inter¬ 
actions  anil  separation  have  been  captured  with  the 
global  PNS  procedure.  Several  differencing  pro¬ 
cedures,  the  stability  of  single  pass  and  global 
relaxation  and  an  analysis  of  the  effective  poisson 
pressure  and  vorticitv  equations  has  been  presented 
in  earlier  studies. 


(  npvrifht  American  InMilute  of  Aeronautic*  and 
Aeronautic*.  Inc..  1983.  All  ri*ht*  reserved. 
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In  the  present  paper,  the  difference  equations 
describing  the  PNS  equations  are  defined  on  a 
partially  staggered  grid.  This  leads  to  a  modified 
interpretation  of  the  discrete  px  approximation. 

A  multi-grid  line  relaxation  procedure  that  reduces 
computer  effort  and  storage  is  outlined.  The 
convergence  properties  and  the  accuracy  of  the 
discrete  system  is  discussed  for  inviscid,  laminar 
and  turbulent  modelling.  The  applicability  of  the 
relaxation  procedure  is  demonstrated  for  inviscid 
flow.  Convergence  studies  are  presented  for  an 
axisvmmetric  boat  tail  configuration.  Incompressible 
and  low  speed  compressible  flow  results  are  then 
discussed  for  laminar  and  turbulent  flow  over  a 
finite  flat  plate  and  a  SACA  0012  airfoil  at  zero 
incidence.  Of  particular  interest  is  the  behavior 
near  the  traiLing  edge,  the  stagnation  point  and 
in  regions  of  separation.  In  addition,  the  impor¬ 
tance  of  the  normal  pressure  gradients  in  these 
regions  will  be  compared  for  laminar  and  turbulent 
flow  cond i t ions. 


2 .  Go- v er  n  i  n g  Eg ua  1 1  on s 

The  governing  PNS  system  considered  here  is 
written  in  a  conformal  coordinate  frame  (£,n)  for 
the  primitive  variables,  pressure  p,  density  r, 
temperature  T,  viscosity  coefficient  p,  and  velo¬ 
cities  u,v.  The  transformed  cartesian  coordinates 
\  =  C(x,v)  aiong  the  surface  and  n  «  n(x,y)  normal 
to  the  surface  are  related  to  the  (x,v)  physical 
coordinates  through  the  transformation  t  *  f(z), 
where  :  =  \+i',  and  z  =  x+iy.  The  metric  h  =  jf*(z)| 
and  hf  -  ih  =  h(f"/f).  The  metric  =  ve  ,  where 

k  =  1  for  axisymmetry. 

Continuity 

(hyu)_  +  (ohvLv)^  =0  (1) 

'.-momentum 

— (;  hv  u\.  +  (fcuvhyf  )^  +  (ouvh)  /h 


As  explained  in  greater  detail  in  earlier 
studies*-  the  version  of  the  PNS  equations  defined 
here  neglects  longitudinal  diffusion  effects  in 
the  ‘"-momentum  equation  (2),  as  well  .is  all  <J  i  t  f  u- 
sion  terms  in  the  normal  (m)  momentum  equation  (5;; 
these  terms,  in  the  appropriate  streamline  coor¬ 
dinate  frame,  are  no  more  important  than  those 
already  neglected  in  (2).  The  continuity  (1)  and 
normal  momentum  equations  are  then  first-order  in  • 
(for  v  and  p),  so  that  two-point  trapezoidal  rule 
differencing  is  applied.  Conventional  three-point 
differences  are  considered  for  -momentum  (2). 
Equations  (1-3)  are  coupled  for  the  variables 
(u,v,p).  The  density  is  obtained  through  an  itera¬ 
tive  process  at  each  \  location,  see  references 
6-8.  This  allows  for  a  simplification  in  the 
treatment  of  the  p  *. «.  an  and  works  quite  well  up  to 
transonic  speeds. 


3.  Difference  Equations  and  Solution  Technique 


The  equations  (L-5)  are  differenced  on  the 
staggered  grid  of  figure  1  as  shown  for  incompres¬ 
sible  flow. 


* 


hv 


hv 


(.  uv)  -  (.  v*"h)^/h  4-  (.  v  ). 


-  -P,  +  —  l^y  .. <bu)  1 
vl  h“ 


(2) 


-momentum 


-p  “  (.  uvhv  )„/hv  +  (.  uvh)./h  -  (.  uv). 


(.  v*"hy  )  -  u2h) 


/h  +  (cu  )  (3) 


hv 


Energy 


u2  + v2 

H  =  C  T  +  -  rr - -  H  =  constant 

p  2 


(4) 


State • V iscos  itv 

p  -  .  KT;  ..  -  m(T) 


O) 


For  the  present  sttidv  only  adiabatic  conditions 
ire  considered. 


Both  first-  and  second-order  accurate  systems  of 


discrete  equations  have  been  developed 


•j-s 


S  i  gn  i  - 


ficantly,  the  first-order  system  given  below,  as 
equations  (6-8),  is,  in  fact,  second-order  accurate 
for  incompressible  inviscid  flow.  This  is  parti¬ 
cularly  important  as  the  PNS  svsten  is  used  to 
evaluate  both  viscous  and  inviscid  regions.  There¬ 
fore,  we  would  expect  that  overall  accuracy  will 
fall  somewhere  between  first  and  second  order. 

This  will  be  demonstrated  in  the  following  sections. 
The  full  second-order  system  is  given  in  references 
3-5,  and  will  not  be  repeated  here.  This  system  is 
more  complex,  requires  many  more  operations  and  is 
verv  sensitive  to  coordinate  t  rails  format  ion , 
regions  of  large  curvature,  or  large  longitudinal 
gradients  as  in  the  trailing  edge  flow.  The 
accuracy  of  tins  system  has  been  tested  lor  a 
semi- inf  ini  to  flat  plate  where  these  effects  are 
negligible.  The  second-order  character  of  the  error 
has  been  confirmed.-*  However,  for  the  reasons 
already  cited,  t  lie  Mf  i  rst  "-order  or  "one  and  one- 
hilf"-«»rder  system*1  ^  is  considered  more  suitable 
for  the  present  applications.  Ihe  equations  are 
shown  here  for  simplicity  in  non-conscrv.it  ion 
cartesian  coordinates  for  tin*  flat  plate  geometry 
and  a  uniform  mesh  in  x  and  v.  In  fact,  non-uniform 
meshes ,  conservation  equations  and  conlnrmal 
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coordinates  are  used  throughout. 


Continuity;  centered  at  0on  figure  1 

((-u)i.j-  (‘u)i-l,j  +  ("u)I,j-l'  <ou)  i-l.J-1  1/2 


+  3?  I(iv)i.j-(^)i.j-i1  =° 


' -momentum:  centered  at  (J) 


(cu)  .  .  (u.  ,-u.  .  .)  +  (i'v)  .  .  zrr~  (u.  .  .-u.  .  ) 

i.j  i.J  i-l.J  i.j  2  Ay  i,]+l  i,j-r 

n-1 

+  p  -  p 

i+i. j  i.j 

.‘.xu .  ,  AxU  .  .  , 

-  - (u.  -u.  .) - K  i  ,> 

fiy2  l.J+l  I.J  Ay2  I.J  1>J~1 


treated  iteratively,  as  in  the  artificial  compres¬ 
sibility  procedure  for  transonic  flow,  it  is 
possible  to  consider  the  fully  elliptic  (.=0)  form 
for  p f .  This  has  been  applied  for  subsonic, ^ 
transonic^  and  supersonic”  cases.  For  supersonic 
flows  with  Af;/y  *  1,  convergence  will  however 

be  slower  than  with  the  exact  treatment  of  (9). 

In  this  procedure  p.  .  and  p.,_  .  are  always 

r  ri , j  r l+l , j 

located  at  the  positions  shown  on  figure  1.  For 

supersonic  regions  the  exact  differencing  would  be 

p,  =  (p?.  ,-pV  .)/A\,  so  that  the  error  is  of  the 


p,  =  (p.  .  . -p.  .  )/A  \ ,  so  that  the  error  is  of  i 

i+1  ,.i  *i,.i 

.  ,  n  n-lx/(1.. 

form  (p.,_  ,-p..,  .  )/A\ ,  i.e.  an  iterative  time 

l+l , J  l+l  ,J 

term  that  vanishes  when  convergence  is  achieved. 

If  the  exact  treatment  (9)  is  used  with 
density  iteration,  then  it  is  appropriate  to 
rewrite  (9) 


Pn"!  .  -  2p*n . 

,  i+L. j _ Lnj 


+  pn 


)  (10) 


n-momentum:  centered  at  (n) 


( c  u )  .  ,  +  (.u). 

i  i.J  t.J-1 


] [v.  ,-v.  ,  ,+v.  .  ,-v.  .  .  .  ] 

i,j  l-l, J  i.j-l  l-l.J-l’ 


where  the  (*)  denotes  the  latest  local  iterative 
value  during  the  nth  global  iteration.  For  w=l, 
at  local  convergence 


(.  v)  .  ,+(nv)  ...  p  .  ,-p .  .  , 

+  ( - — - - laJ~l)/v  _v  - LrX_l 

1  2 Av  MVi,j  Vi,j-lJ+  Ay 


The  equations  (6-8)  represent  a  coupled  system  for 


the  vector  (v  J  .  As  seen  from  the  unknown 
vp  '  i.j 

bracketed  vector  in  figure  1,  on  the  staggered  grid 
the  pressure  p.  .  is  evaluated  one  grid  point 


upstream  of  the  unknown  velocity  u ^  .  This  re¬ 

flects  the  upstream  influence  or  elliptic  effect 
of  the  pr  term.  The  superscript  n-1  for  P^+,  j  in 

(7)  denotes  the  value  of  the  previous  global  itera¬ 
tion  or  sweep.  The  current  iteration  is  designated 
by  n;  although,  this  superscript  has  been  dropped 
throughout.  It  is  significant  that  only  the 

pressure  term  p"+j  .  is  relaxed  during  the  itera¬ 
tion  procedure.  Therefore,  only  p.  .  need  be 

i.j 

stored.  For  separated  flow  the  convective  terms 

are  upwinded  so  that  u!  ,!  .  and  v"^  ,  are  also 

i+1  .J  i+1 . J 

required,  but  only  in  the  limited  regions  of 
reverse  flow. 

The  treatment  of  p^  in  (7)  as  (P|+|  j — P j*  j )  / A r 

corresponds  to  a  fully  elliptic  approximation. 

From  previous  studies  by  the  authors,  as  well  as 
the  eigenvalue  analysis  of  reference  13,  this 
description  is  exact  only  for  Incompressible  (w=0) 
flow.  For  subsonic  flow,  the  structure  discussed 
in  the  Introduction,  where  p,  is  of  the  form 


pT.  ’  •fpi.j"pi-lIj>/‘'5  +  (1_“)(piIl,fpi,))'Ar 

( 

should  apply  if  the  density  is  coupled  with  the 
(u,v,p)  calculation.  However,  if  the  density  Is 


In  this  procedure,  the  pressure  at  the  u.  .  location 
of  figure  1  becomes  wpj  j  +  (l-..)p^+^  and  global 

convergence  is  accelerated  for  supersonic  points. 
Both  methods  have  been  applied  successfully  here 
and  in  references  6-8.  Equation  (!0)  is  direct  1 v 
obtainable  from  (9),  if  P^+j  •  O  approximated  by 


n 

Pi-l,j 


so  that  the  two  approaches  differ 


only  in  truncation  error. 

The  density,  temperature  and  viscosity  coef¬ 
ficient  are  obtained  from  (4)  and  (5).  A  linear 
u(T)  relationship  has  been  assumed  herein.  When 

f"\ 

these  quantities  have  been  updated  the  vector  (  v  I 

V  P  '  i 

is  re-evaluated  until  some  specific  degree  of 
convergence,  usually  10“*  in  sequential  steps,  is 
achieved.  The  nonlinear  svstem  is  quas  i  - 1  i  near  i  cd 
in  a  standard  manner and  nonlinear  terms  are 
updated  at  the  same  time  as  the  density,  'iscositv, 
etc.  After  final  convergence  at  location  i,  the 
procedure  continues  to  i+1,  etc. 

The  tridiagonal  linearized  system  (h-8)  i 

solved  for  (v  J  bv  the  standard  1. 1’  detemp.-i- 

'  P  '  i.j 

tion.  The  total  line  relaxation  procedure  rent  inn*  ' 
to  the  right  hand  boundary.  The  (n-1)  terms  art- 
then  updated  and  the  marching  process  is  repeated. 
The  solution  procedure  is  terminated  when  tin 
changes  in  maximum  pressure  and  skin  fri>  t ion 
between  global  iterations  is  less  than  1  n“'4 .  I  in* 
relaxation  is  considered  most  efficient  tor  the 
present  analvsis  as  it  minimizes  storage  requir* 
ments.  More  efficient  procedures,  such  as  AIM  i 
CSIP  would  require  increased  storage  tor  the 
velocities  as  well  as  the  pressure.  Since  *-ur  c>ii 
is  to  reduce  storage  requirements  ter  future  thi*«  - 
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M 


-  ,\ 

*  \ 


& 


* 


ti: 


K>.-. 


dimensional  applications,  simple  line  relaxation 
with  a  multi-grid  corrector,  to  he  discussed 
shortlv,  was  chosen. 


Boundary  Conditions 


For  the  primitive  variable  system  (b-8)  on  the 
staggered  grid  of  figure  1,  the  following  boundary 
conditions  apply: 


At  the  inflow 


O' 


v . ( \q, i)  =  V(n) . 


V(-) 


l  ( -  )  and 
(■  )  +  C’  (•  ) ,  where  .  (■  ) 


u(-.0,O 


is  the  inflow  vorticitv;  for  uniform  conditions 
l1  I  and  V  "  0  (zero  vorticitv).  A  condition  for 
the  pressure  is  not  required  for  incormross ihle 
flow.  As  seen  from  figure  1,  the  inflow  pressure 
i s  c  a 1 culated  during  the  first  marc  1 1  ing  step . 


At  the  upp e r  surf a c e 


M 


u  =  1, 


1  ; 


i.t*.,  free-stream  conditions  are  applied.  This  re¬ 
quires  that  • w  be  sufficient  tv  large,  e.g.,  outside 
M  '  3-5 


of  the  domain  of  the  triple  deck  interaction. 

The  outer  boundary  could  he  moved  (‘loser  to  the 
surface  bv  applying  conditions  similar  to  those  of 
interacting  boundary  layer  theory;  however,  for  the 
present  calculations,  the  fully  inviscid  evaluation 
was  included.  A  boundary  condition  on  v  is  not 
requ i red. 


At  the  outflow 


PC 


)  or  derivative  p 


Onl_y  t he_  pressur e 
j  , »  )  a re  pr escj  i be d . 


For 


the  present  calculations,  both  p  =  0  and  p.  =0 
were  applicable.  There  are  only  slight 
differences  in  the  solutions. 


•v 

y< 

At  tin 
/')  =  0 

wall  •  =  0:  For  viscous  flow  u(".,0)  = 
is  s|us  i  f  it-ti.  For  inviscid  flow  only 

* 

V  (  ' 

,io  -  o 

is  required.  A  boundarv  condition  on 

the 

pressure  is  not  required. 
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imntlv  pressure  boundarv  conditions  an 
ied  at  the  upper  and  outflow  boundaries. 
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in  reference-.  *,5,  the  equations  (b-8) 
inipul.ited  t  »  provide  an  effective  poisson 
•»  lver;  the  remaining  pressure  boundarv 
m,  Neirumi  tvpe,  are  implicitly  imposed 
in!,  r  red  dire,  t  1  v  from  the  original 
.  e-,.i.  i  f  i  m  !  '• .  [here  is  no  need  to  npplv 

•r*  n*  !  il  !  >f  the  oisson  equation  and 

indiredlv,  as  is  commonlv 
>«  ompre*  s  i  bi  e  flow  prore- 
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1  the  rate  of 
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lurUe.Jiv.  It  is  interest- 
, .  1  ,  wh i  -  h  1 eads  to 


t-*r  ■.iiii'li-  ".weep  methods, *" 


is  a  .  .'iiviTiyivc  factor  lor  lb 
procedure.  As  noted  in  the  previa 
from  the  eigenvalue  anal  vs  is  ^  *  for 
it  can  be  inferred  that  the  nmvvri1 
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will  be  of  the  form  J  —  (1-  'll-A(-) 
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implies  mi, re  rapid  convergence  with  large 
of  high  subsonic  or  supersonic  Mach  uumbi 
■since  ..  ■  1  as  M  *  I. 
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In  order  to  improve  the  convergence  rate,  a 
one-dimensional  (in  :)  multi-grid  procedure 
following  the  full  approximation  scheme  of 
reference  (14)  has  been  applied.  This  approach 
is  less  efficient,  in  terms  of  storage  requirements, 
than  a  full  two-dimensional  mult; -grid  application; 
however,  relaxation  for  p  is  only  in  the  '-direction 
and  the  calculation  is  fully  implicit  in  •.  Also, 
a  highly  non-uniform  mesh  is  required  to  accurately 
describe  the  boundary  layer,  triple  deek  and 
inviscid  regions.  Hi  is  is  not  ideal lv  suited  for 
interpolation  required  bv  the  multi-grid  procedure. 
Of  particular  note,  error  transfer  from  coarse  to 
fine  grids  is  applied  only  for  the  pressure,  even 
in  regions  of  reversed  flow  where  the  velocities 
n-1  ,  n-1 


•i+i.j  and  *i+i, j 


are  also  relaxed.  These  values 


are  fixed  on  each  grid  from  the  previous  global 
(n-1)  iteration.  Since  the  separation  zones  are 
of  considerably  different  extent  on  each  grid, 
the  transfer  process  for  the  velocities  was  con¬ 
sidered  to  be  a  source  of  possible  difficulty.  The 
present  procedure  works  quite  well  for  most  of  the 
problems  considered  here  and  is  discussed  in 
greater  detail  in  references  4,5. 


The  effects  of  the  multi-grid  procedure  on  the 
convergence  properties  are  shown  in  figure  2  for 
a  trough  geometry  discussed  previously  in 
references  4,5.  Similar  results  were  found  for 
the  flat  plate  configuration.  When  large  regions 
of  separated  flow  occur,  the  coarse  grid  corrections 
are  less  effective  and  the  multi -grid  procedure 
appears  to  enter  a  limit  ryi  le.  For  this  reason, 
the  coarse  grid  solutions  were  simply  applied  as 


v  U1  I  MUM  I  MU!  I  1 

\ 


mi ii  1 1  (>'ip  ( i  ppup  nn  i  ini  m  u.'iii 

■  I  I  I  II  I  I  I  I  I  I  l  |  1  l  :  ■  |  II 


II! 


!  II' 


MO  MI  I,t  PMAl  l  il  PA  I  1 1  m. 


Fig.  2.  Convergence  Kate  with  Multi -Or  id  Scheme 
(Trough  •  =  -0.015,  Re =  8 x  104), 
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fine  grid  initializations.  This  uni-directional 
process  was  still  a  considerable  improvement  over 
direct  calculations  on  the  finest  grid.  Figure  2 
depicts  the  significant  reduction  in  iteration 
count  with  the  multi-grid  procedure.  Figure  3 
depicts  the  sensitivity  of  the  pressure  solution 
to  mesh  width  for  the  trailing  edge  triple  deck 
interaction.  The  skin  friction  is  much  less  grid 
dependent.  This  behavior  was  found  for  other 
geometries  as  wellJ*®  Finally,  the  convergence 
rate,  for  a  fixed  >,  mesh,  as  a  function  of  the 
outer  boundary  location  y^  is  shown  in  figure  4. 
The  results  are  consistent  with  the  prediction 
that  convergence  will  deteriorate  as  (&r/y  ) 
decreases.  With  the  multi-grid  correction "however 
(figure  2),  this  sensitivity  is  greatly  reduced 
and  the  parameter  (A%/y  )  is  less  critical. 
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Fig.  4.  Convergence  Kate  for  Different  (.'.x/y„ ) 
(Semi-Infinite  Flat  Plate). 


6.  Solutions 


(i)  Inviscid  Flow  -  Boattail 


I  Re  -  105 


i  C0.A8SE  GRID  fO . 1 > 


*  FINE  GRID  (0.0U 


Fig.  3a.  Flat  Plate  Trailing  Edge  Pressure  on 
Two  Grid  Levels  (Laminar). 


Since  the  PNS  equations  are  applied  here  for 
the  evaluation  of  the  complete  viscous/invisc id 
interaction,  the  effectiveness  of  the  relaxation 
procedure,  for  the  discrete  system  (1-5),  for  in¬ 
viscid  flow  is  important.  For  supersonic  flow  with 
a  fitted  shock  boundary,®  the  inviscid  region  is 
relatively  small;  however,  for  subsonic®  or  tran¬ 
sonic^  flows,  these  regions  can  be  quite  extensive 
and  the  overall  convergence  rate  of  the  PNS  proce¬ 
dure  will  be  dependent  on  the  "inviscid"  convergence 
properties.  In  order  to  test  the  present  method, 
solutions  were  obtained  for  the  inviscid  flow  (slip 
boundary  conditions)  over  an  axisymmetric  boattail 
configuration.  The  results  are  shown  in  figure  5, 

without  a  multi-grid  correction,  and  compared  with 
potential  flow  solutions  obtained  with  the  CSIP 
procedure. The  initial  conditions  are  somewhat 


COARSE  GRIG  (0.1) 


a  -  POTENTIAL  SOLUTION 
»  -  GLOBAL  PNS 


»  FIXE  GRID  (0.01)  -= 


:  .o-2- 

x 


m 


E 

R  -3 
R  t0  3- 
0 
R 


-r  VTTVTTf - V  }  TI-T 

t  e 

NO.  OF  GLOBAL  ITERATIONS 


r  t~!  rni 
1 1 


RATE  OF  CONVERGENCE  (BOATTAIL  BF.TAC  -  7.891  DEG  INVTSCTD5 


Fig.  3b.  Flat  Plate  Trailing  Edge  Skin  Friction 
on  TVo  Grid  Levels  (Laminar). 


Fig.  5.  Convergence  Rate  for  Inviscid  Flow  On 
Boattail  (Juncture  Angle  15.8°). 
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different  and  therefore  the  potential  solutions 
start  with  a  slightly  larger  error.  The  figure 
clearly  shows  that  the  convergence  rate  of  the 
line  relaxation  procedure  for  the  pressure  des¬ 
cribed  herein  is  not  significantly  different  than 
that  for  the  potential  equation  and  the  CSIP.  This 
is  an  important  and  somewhat  surprising  result. 
Moreover  as  seen  in  figure  6a,  the  solution  is 
essentially  second-order  accurate.  This  was  dis¬ 
cussed  previously;  further  analysis  is  given  in 
references  4,5  where  it  was  shown  that  on  the 
staggered  grid  of  figure  1,  the  inviscid  discrete 
equations  for  the  pressure  and  velocities  are  in 
2  2 

error  by  0(Ar  ,  An  ),  when  the  vorticity  is 
defined  by 


v. -  v.  +  v, .  ,  -  v.  .  , 

_  1+1  »j _ i,j  i+i,  j-i  i,j-i 

i,.i  2t.:, 


u.  .  -  u.  .  . 
_  l.J  1.1-1 


These  results  confirm  the  validity  and  accuracy 
of  the  relaxation  procedure  for  inviscid  regions. 
An  accuracy  plot  for  the  viscous  boat  tail  dis¬ 
cussed  previously is  given  in  figun  6b.  Note 
that  the  exponent  1.8  is  obtained.  Similar  ^ 

results  were  found  for  the  trough  configuration. 
This  indicates  that  the  overall  accuracy  of  the 
system  (6-8)  is, as  predicted,  somewhere  between 
first  and  second  order. 

(11)  Finite  Flat  Plate  -  P N S 


Solutions  for  the  full  PNS  system  are  given 
in  references  3-5  for  the  laminar  incompressible 
flow  over  a  finite  flat  plate.  The  pressure 
coefficient  given  in  figure  3  for  two  grids  is 
shown  in  figure  7.  These  results  are  repeated 
here  in  order  to  demonstrate  the  excellent  agree¬ 
ment  obtained  with  interacting  boundary  layer 
solutions  for  this  geometry.  The  trailing  edge 
pressure  and  skin  friction  are  predicted  vc-rv 
accurately  and  for  the  staggered  grid  the  minimum 
pressure  is  obtained  at  the  trailing  edge  of  the 
plate.  Due  to  the  discontinuity  in  pressure 
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Fig.  6a.  Accuracy  Plot  for  Inviscid  Flow  on 
Boattail  (Juncture  Angle  15.8°). 
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Fig.  6b.  Accuracy  Plot  for  Boattail  (.Juncture 
Angle  30°)  Laminar  Flow  Re  =  1000. 


Fig.  7a.  Flat  Plate  Trailing  Edge  Pressure 
(Laminar  Flow). 
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Fig.  7b.  Flat  Plate  Trailing  Fdge  Skin  Friction 
and  Wake  Centerline  Velocity  (Laminar). 


Vs-: 


►  V 


1 


199 


gradient  at  the  trailing  edge,  the  accuracy  factor 
was  1.35,  a  value  somewhat  less  than  that  found 
for  the  boattail  and  trough. 

For  the  present  study  we  are  particularly 
concerned  with  the  trailing  edge  behavior  for  > 

turbulent  flow  conditions.  Th>-  two-laver  Cebeci-  ' 

1  f  N 

Smith  eddy  viscosity  model10  has  been  assumed  for 

the  Reynolds  stress  term  in  the  '-momentum  equa-  * 

■  |) 

tions  on  the  plate,  and  a  modified  version  given  ( 

in  reference  16  has  been  applied  in  the  wake. 

Reynolds  stresses  have  been  neglected  in  the  ; 

n-momentum  equation.  The  two-laver  model  should  o 

be  adequate  for  this  unseparated  flow  problem.  N 

Limited  interacting  boundary  layer  results  have 
been  presented  in  reference  17.  A  definitive 
triple  deck  theory  has  not  been  completed  for  this 
flow.  One  of  the  open  questions  concerns  the 
influence  of  normal  pressure  gradients. 

Solutions  for  the  pressure  coefficient,  skin 
friction  and  wake  centerline  velocity  are  shown  in 
figures  8a  and  8b.  The  centerline  velocity  and 
pressure  results  are  in  reasonable  agreement  with  in¬ 
teracting  boundary  layer  and  experimental  vrlues. 

The  normal  pressure  gradients  at  the  trailing  edge 
(figures  9a  and  9b),  however,  are  quite  different 
than  those  found  for  laminar  flow,  (figure  9a). 

In  the  latter  case  virtually  all  of  the  pressure 
change  occurs  outside  of  the  boundary  layer,  i.e., 
in  the  inviscid  region  or  outer  deck  of  the  triple 
deck  structure.  For  the  turbulent  flow,  approxi¬ 
mately  60"i  of  the  pressure  change  occurs  inside 
the  boundary  layer  (figure  9b).  Therefore,  the 
longitudinal  pressure  distribution  (figure  8a)  at 
the  surface  (centerline)  is  considerably  different 
than  that  at  the  boundary  layer  edge.  This  clearly 
has  some  implications  with  regard  to  the  applica¬ 
bility  of  conventional  interacting  boundary  layer 
methods  that  neglect  these  gradients.  Similar  ^  ^ 
effects  have  been  found  for  the  boattail  geometry.  ’ 

It  is  surprising  that  the  values  of  the  centerline 
velocity  (figure  Cb)  are  in  such  good  agreement. 
Further  study  on  the  effect  of  normal  pressure 
gradients  for  turbulent  flow  is  certainly  suggested. 
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Fig.  9a.  Normal  Pressure  Distribution  Near  Trailing 
Edge  of  Flat  Plate  (Laminar,  Re  =  10^). 
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Fig.  8a.  Flat  Plate  Trailing  Edge  Pressure 
(Turbulent,  Re  *  6.5  x  10^). 


Fig.  9b.  Normal  Pressure  Distribution  Near  Trailing  Edge 
of  Flat  Plate  (Turbulent,  Re  =  b.  r)  x  10^). 


(Hi)  NACA  0012  Airfoil  -  PNS 

Laminar  and  turbulent  flow  over  a  NACA  0012 
at  zero  incidence  is  considered  for  incompressible 
and  subsonic  (M  *0.3)  flow.  The  two-layer  eddy 
viscosity  model  is  applied  once  again.  This  should 
be  suitable  for  non-separated  flow;  however,  for 
the  separated  flow  solutions,  the  adequacy  of  this 
model  could  be  questioned.  For  the  mild  separation 
regions  occurring  with  the  airfoil  geometrv,  the 
errors  incurred  with  this  closure  model  should  be 
relatively  small.  This  is  not  true  for  regions 
with  massive  separation  zones.  The  metric  h  in 
the  governing  equations  (1-3)  is  obtained  from  a 
conformal  mapping  procedure  described  in  previous 
studies. The  grid  is  shown  in  figure  10. 
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Fig.  111.  NACA  0012  Airfoil,  Laminar  w 
(Re  =  12,500). 
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Fig.  10.  NACA  0012  Airfoil  Grid  Close  to  the  Body 
(Physical  Plane). 


Laminar  flow  solutions  have  been  obtained  for 
Reynolds  numbers,  Re  =  5000,  10000,  12500.  For 
the  larger  value  of  Re  a  small  separation  region 
appears  near  the  trailing  edge,  figures  11a,  11b. 

A  detailed  study  of  the  front  stagnation  point  for 
a  series  of  successively  finer  grids  has  provided 
evidence  that  the  PNS  approximation,  in  an  appro¬ 
priate  body  fitted  coordinate  system,  can  accurately 
evaluate  the  flow  in  this  region.  Recall,  that 
the  exact  Navier-Stokes  equations  for  stagnation 
point  flow  do  not  include  diffusion  effects  along 
the  surface  or  diffusion  terms  in  the  normal 
momentum  equation.  Comparisons  of  the  present 
PNS  stagnation  point  shear  stress  with  the  exact 
Navier-Stokes  value  of  1.25  are  given  for  several 
grids  in  figure  12.  For  the  finest  two  grids,  the 
agreement  is  quite  good. 

Turbulent  flow  solutions  are  shown  in  figures 
13a  and  13c* for  Re  -  5.35 x  10 3.  Based  on  experi¬ 
mental  data,^®  transition  has  been  prescribed  .it 
a  distance  along  the  airfoil  x/c  =  0.4,  where  c  is 
the  chord  length.  Separation  does  not  occur.  The 
agreement  with  the  data  is  reasonable  (figure  13a) 
considering  the  fact  that  the  grid  resolution  is 
clearly  inadequate  for  the  region  x/c  0.4,  where 
the  flow  is  laminar  at  the  same  very  large  Reynolds 
number.  Finally,  the  normal  pressure  gradient  for 
the  laminar  and  turbulent  flows  is  shown  in 
figure  1 3b, C. Behavior  similar  to  that  found  for 
the  trailing  edge  is  obtained  for  the  airfoil 
conf  igura t ion. 
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Fig.  13a.  NACA  0012  Airfoil  Turbulent  Flow 
(Re  =  5. 35  x  105) . 


Fig.  13b.  Normal  Pressure  Distribution  Near  Trailing  Edg 
of  NACA  0012  Airfoil  (Laminar,  Re  =  10^). 
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Fig.  13c.  Normal  Pressure  Distribution  Near  Trailing  Edge  of 
NACA  0012  Airfoil  (Turbulent,  Re =  5.3SX  10^). 


(iv)  Other  Geometries 

Solutions  with  large  separation  regions  have 
been  presented  for  the  boattail  geometry.-*  Both 
laminar  and  turbulent  flows  have  been  considered. 
Finally,  for  subsonic  and  transonic  flow,  laminar 
solutions  over  a  cone-cvl inder-boa t ta i 1  configura¬ 
tion,  for  which  separation  occurs,  have  been  pre¬ 
sented  in  references  6,  7.  Initial  results  for 
supersonic  flow  over  a  cone  have  been  given  in 
reference  8.  Additional  supersonic  studies  will 
be  presented  in  future  papers. 


Laminar  and  turbulent  solutions  for  finite 
flat  plate  and  NACA  0012  airfoil  geometries  have 
been  obtained  with  a  multi-sweep  PNS  relaxation 


procedure.  The  method  is  unconditionally  stable, 
departure  free,  separation  singularity  free  and 
describes  separation,  stagnation  point  and  trailing 
edge  behavior.  With  the  difference  equation  given 
on  a  partially  staggered  grid,  the  inviscid  equa¬ 
tions  are  second-order  accurate  and  the  full  PNS 
system  is  somewhere  between  first  and  second-order. 
The  convergence  rate  for  inviscid  flow  is  comparable 
to  that  obtained  for  the  full  potential  equation. 

For  fine  meshes  or  large  values  of  the  upper 
boundary  y^,  the  convergence  rate  slows  considerably 
A  one-dimensional  multi-grid  procedure  for  the 
pressure  has  been  applied  in  order  to  alleviate 
this  difficulty.  The  solutions  are  in  good  agree¬ 
ment  with  previous  results  or  data  and  the  effects 
of  normal  pressure  gradients  for  turbulent  flows 
are  shown  to  be  significant  even  within  the  boundary 
laver.  The  method  is  applicable  for  flows  with 
strong  pressure  interaction  and  large  separation 


regions  when  the  difference  equations  are  written 
in  an  appropriate  streamline  coordinate  system. 
Other  studies  for  cone,  trough,  boat  tail,  cone- 
cvl inder-boattail  and  base  flows  have  established 
the  procedure  as  applicable  for  subsonic, 
transonic  and  supersonic  flow. 
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